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^ ^ Abstract 

Q^' We identify reactions which destroy ^Be and ^Li during big bang nucleosynthesis (BBN) in the 

scenario of BBN catalyzed by a long-hved sub-strongly interacting massive particle (sub-SIMP or 



o 



X particle). The destruction associated with non radiative X captures of the nuclei can be realized 



■ only if the interaction strength between an X particle and a nucleon is properly weaker than that 

> ; 

. between two nucleons to a degree depending on the mass of X. Binding energies of nuclei to an 



X particle are estimated taking the mass and the interaction strength to nuclei of the X as input 
CN ' parameters. Nuclear reaction rates associated with the X are estimated naively, and adopted in 

calculating evolutions of nuclear abundances. We suggest that the ^Li problem, which might be 
associated with as-yet-unrecognized particle processes operating during BBN, can be solved if the 
^ ' X particle interacts with nuclei strongly enough to drive ^Be destruction but not strongly enough 



to form a bound state with ^He of relative angular momentum L = 1. Justifications of this scenario 
by rigorous calculations of reaction rates using quantum mechanical many-body models are highly 
desirable since this result involves many significant uncertainties. 
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I. INTRODUCTION 



The standard big bang nucleosynthesis (BBN) model predicts primordial light element 
abundances which are more or less consistent with abundances inferred from observations 
of old distant astronomical objects. Deviations from the standard BBN (SBBN) model are, 
therefore, constrained if predicted abundances in theoretical models change from those in 



the SBBN. Constraints on the existence of 
nuclei by strong force [l - 4 1 or Coulomb force 



■ong- 



ived exotic particles which interact with 



22| have been derived as well as those on the 



decay of long-lived exotic particles into standard model particles which have electromagnetic 



or hadronic interactions 



23 



50|. 



A prominent problem relating to the abundances predicted in the SBBN model and in 



ferred from observations is lithium problem 
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52| . Primordial lithium abundances are 
inferred from measurements in metal-poor halo stars (MPHSs). Observed abundances 



are roughly constant clS db function of metallicity 
The theoretical prediction 
^Li/H=(5.24i°:^i) X 10-10 
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-l57| at ^Li/H= (1 - 2) x 10"^°. 
3y__the SBBN model is, however, a factor of 2-4 higher, e.g., 
when its only parameter, the baryon-to-photon ratio, is 



deduced from the observation with Wilkinson Microwave Anisotropy Probe (WMAP) of 
the cosmic microwave background (CMB) radiation 59|. This discrepancy indicates some 
mechanism of ''Li reduction having operated in some epoch from the BBN to this day. One 
possible astrophysical process to reduce ^Li abundances in stellar surfaces is the gravitational 
settling in the model including a combination of the atomic and turbulent diffusion 60, l61 |. 



The precise trend of Li abundance as a function of effective temperature of stars in the 



62|. 



metal-poor globular cluster NGC 6397 is, however, not reproduced theoretically 

^Li/^Li isotopic ratios of MPHSs have also been measured spectroscopically. The ^Li 
abundance as high as ^Li/H~ 6 x 10 "^^ was suggested 52], which is about 1000 times 



higher than the SBBN prediction 



1271]. Convective motions in the atmospheres of MPHSs 



could cause systematic asymmetries in the observed line profiles and mimic the presence of 
^Li 63[. A few or several MPHSs, however, have high ^Li abundances larger than levels 
caused by this effect 6J]. This high ^Li abundance is a problem since the standard Galactic 
cosmic ray nucleosynthesis models predict negligible amounts of ^Li yields compared to 
the observ ed le vel in the epoch corresponding to the metallicity of the stars, i.e., [Fe/H] 
< -2 y |l28]. 



2 



The possibility that the ^Li and ^Li problems stem from uncertainties in nuclear reactions 



66| un 



67 



ess there remain to be observed 
The ^Li reduction needs a 
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used in theoretical BBN calculation is unlikely 
new resonant states contributing to ''Be destruction 
destruction mechanism of ^Be during or after BBN and before stellar activities since the ^Li 
nuclei observed in MPHSs are thought to have originated from the electron capture process 
of ^Be which is produced in the BBN. 

Some particle models beyond the standard model include heavy (m ^ 1 GeV) long-lived 
colored particles. The scenarios, i.e., split supersymmetry jgQ, 70|, weak scale supersymme- 



try with a long-lived gluino 
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- |73| or squark [7^ as the next-to-lightest 
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supersymmetric 



76|, may be tested 



particles, and extended theories with new kinds of colored particles 
in experiments such as the Large Hadron Collider. The heavy colored particles would be 
confined at temperatures below the deconfinement temperature Tc ~ 180 MeV inside exotic 
heavy hadrons, i.e., strongly interacting massive particles (SIMPs) which we call X parti- 
cles [77|. Their thermal relic abundances after the freeze-out of annihilations depend on the 
annihilation cross sections, and theoretical estimates predict various values which extend 
over more than several order of magnitude at the heavy mass limit 781 ]. 

If the annihilation cross section is not different from a typical value for strong interaction, 
i.e., cr ~ 0{GeV~^y, however, the final abundance of X particles can be derived under the 
assumption that their abundances are fixed when the annihilation rate becomes smaller than 
the Hubble expansion rate of the universe 77|]. The relic abundance can then be written 

Nx 



^ 15 g'J^ mV^ 
^ 9*s crTo^^mpi 



(1) 



where Nx is the number density of the X particle, s = 27r^(7*sT^/45 is the entropy density 



with g 



10 the total number of effective massless degrees of freedom in terms of en- 



tropy [79[ just below the QCD phase transition, g^, is the total number of effective massless 
degrees of freedom in terms of number [t^, m is the mass (m ^ 1 GeV) of the heavy 
long-lived colored particles, a is the annihilation cross section of the X particle, is the 
temperature of the universe at which the X-particles are formed, and mpi is the Planck 
mass. The number abundance of the Xs with respect to that of baryons is then 

iV„ _/ n \ 1/2 / rr, \ 1/2 / \ ^3/2 



0.5 X 10' 



9* 



1/2 / \l/2 

10.757 VTW 



180MeV/ 



a 



(2) 



where is the number density of baryons, and ~ 140 MeV is the mass of pion. The 



thermal relic abundance is inversely proportional to the annihilation cross section which 
depends on the particle theory. In addition, there might be a nonthermal production of 
long-lived colored particles which is not directly related to the thermal production. The 
final abundance of the X is, therefore, uncertain. So we consider the X abundance as a free 
parameter in this paper. 



Observational constraints on hypothetical SIMPs have been studied 80|-|82|. Effects of 
exotic neutral stable hadrons on BBN were studied in Ref. jl|. The authors assumed that the 
strong force between a nucleon and a exotic hadron (X) is similar to that between a nucleon 
N and a A hyperon. In addition, new hadrons are assumed to be captured in a bound state 
of ^He plus X after BBN. Based upon an analytic estimation, they suggested that beryllium 
has the largest number fraction Ax/ A of bound states with the hadrons among the light 
elements produced in BBN, where the A and Ax are a nuclide A and a bound state of A 
with a hadron X. Mohapatra and Teplitz jsl estimated the cross section of X capture by 



4 



He, and suggested that a large fraction of free X particles would not become bound into 



light nuclides and remain free contrary to the previous suggestion p| . In deriving the result 
of those two studies, it has been assumed that exotic hadrons interact with normal nuclei 
by typical strengths of strong interaction and implicitly assumed that its mass is about that 
of A hyperon, i.e., mx ~ 1.116 GeV 



Effects on BBN of long-lived exotic hadrons of m ^ 1 GeV have been studied recently j4|. 
The authors have assumed that the interaction strength between an X particle and a nucleon 
is similar to that between nucleons. Rates of many reactions associated with the X particle 
were estimated, and a network calculation of the nucleosynthesis including effects of the X 
was performed. The constraint on the decay lifetime of such X particles, i.e., tx ^ 200 s 
was derived from a comparison of calculated abundances with observational abundance 
constraints of light elements. 

Two interesting predictions of the model jj] is signatures of the X particles on primordial 
abundances which should be seen in future astronomical observations: 1) ^Be and B can 
be produced in amounts more than predicted in the SBBN. Future observations of Be and 
B abundances in MPHSs may show primordial constant values originating from the BBN 
catalyzed by the X particle. 2) The isotopic ratio ^°B/^^B tends to be very high. This is 



different from predictions of other models for boron production, i.e., the cosmic ray nuc 



osynthesis (^''B/^^B~ 0.4 84l-l86|) or the supernova neutrino process (^°B/^^B<^ 1 871. l88|). 



e- 
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They concluded that the ^Li or ^Li problems is not solved under their assumption. 

Since interactions between long-lived exotic hadrons X and a nucleon are not known 
as well as their masses, we are investigating effects of such particles in various cases of 
interaction strengths and masses. We found on the way a new possibility that reactions 
associated with the X particle reduce ^Be abundance and that the ^Li problem is solved. 

In this paper we report details of the destruction mechanism of ^Be in the presence of 
the X particle. We carry out a network calculation of BBN in the presence of a long-lived 
SIMP X^ of a zero charge taking the mass and the strength of interaction with a nucleon as 
characterizing parameters. In Sec. [Ill assumptions on the X^ particle, estimations for binding 
energies between nuclei and an and rates of important reactions are described. Effects 
of the X^ decay inside X-nuclei are not considered in our model. They should be addressed 
in the future. In Sec. Illll the destruction processes of ^Be and ''Li are identified. With results 
of the network calculations of BBN, we delineate the parameter region in which the ^Be and 
^Li destructions possibly operate. If the X^ particle interacts with nuclei strongly enough 
to drive ^Be destruction but not strongly enough to form a bound state with ^He of relative 
angular momentum L = 1, then it might solve the ^Li problem of standard BBN. In Sec IIVI 
conclusions of this work are summarized, and this model for ^Li reduction is compared with 
other models. 



II. MODEL 



A strongly-interacting massive particle (SIMP) X of spin zero and charge zero is assumed 
to exist during the epoch of BBN. Its mass is one parameter since it is not known a priori 
at the moment. Two types of nuclear potentials between an X^ and a nucleon (-^A^) are 
considered in this study. One is the Gaussian type given by 



v{r) = vq5 ex.]) — (r/rg) 



(3) 



1291], and the interaction strength is varied 



where Vq = —72.15 MeV and tq = 1.484 fm 
by changing 6, the second parameter. When the 6 equals unity then the binding energy of 
deuteron, i.e., 2.224 MeV is obtained. 



The potential between an X^ and a nuclide A {XA) is given by 

V{r) = f v{x)p{r') dr\ 



(4) 



where r is the radius from an X to the center of mass of A, r' is the distance between the 
center of mass of A and a nucleon inside the nuchde A, x = r + r' is the distance between 
the X and the nucleon, and p{r') is the nucleon density of the nucleus which is generally 
distorted by potential of an X^ from the density of normal nucleus. Under the assumption 
of spherical symmetry in nucleon density, i.e., p(r), the potential is written in the form of 



V{r) = ttvqS- 



r Jo 



dr'r' p{r' 

J\2 



X |exp 

Another potential is a well given by 



^ _ ^. 





(r + r')2' 




— exp 




I 









v^{r) = wow'^w (for < r < Vq^), 



(5) 



(6) 



and fwl'^ow < r) = 0. Parameters are fixed to be wqw = —20.06 MeV and tqw = 2.5 fm. 
In order to make a comparison between the two potential cases easy, integrals, i.e., I = 
J v{r)dr for both cases are made equal when 6 = 6^^ = 1. This integral is a characteristic 
quantity which is related to binding energies. The requirement of equal integral values and 



an assumption of row = 2.5 fm 



130| leads to 



^^0 



4 



TT / ro 
^0^ 



— = -20.06 MeV. 



(7) 



The XA potential is given by 

V{r)^ = 2'kvq^5^- I dr'r'p{r')\ dxxH{x-ro„), (8) 

r JO J|r--r'j 

where H{x) is the Heaviside step function. 

As a crude assumption, the nucleon density p(r) is approximately given by the undistorted 
one for normal nucleus. The folding technique to derive XA potential from the XN potential 
[equation (jl])] then does not exactly yield the XA potential in any mathematically-rigorous 
method of calculation. In order to derive precise results of nuclear structures or energy 
levels, all nucleons as well as an X^ and all interactions among them need to be taken 
into account with many-body quantum mechanical calculations. Since such calculations are 
unrealistically difficult, three or four-body models for an X^ particle and nuclear c 



usters 

composing the nucleus should be utilized as were done in the case of hypernuclei [89[. The 
assumption taken here only provides some reasonable estimate of what an effective X- 
nucleus reaction might look like with all the nuclear degrees of freedom frozen out. The 
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folding procedure, however, might produce a more useful approximation than in the nuclear 
case since the does not participate in the Pauli principle among nucleons. 

The nucleon density of nuclei with mass number A>2 is assumed to be Gaussian, i.e, 

p(r) = p(0)exp[-(r/6)2], (9) 

where p{0) = A'K~^f'^h~^ is the nucleon density at r = and satisfies the normalization 
/ p{r)dr = A, with A the mass number. The parameter for the width of density, i.e., b, is 
related to the root mean square (RMS) nuclear matter radius which should be determined 
from experiments, i.e., b =^J2/3r^^. 

The XA potential in the case of the Gaussian XN potential, i.e., equation (|5]), is then 
simply written as 

A. Nuclear Binding Energies 

The BBN catalyzed by the X particle is significantly sensitive to binding energies of 
nuclei to an X^ particle (X- nuclei). The binding to X particles changes the relative energies 
of initial and final states, and may even change the sign of the Q-value. j4|. Binding energies 
and eigenstate wave functions of X-nuclei are computed taking into account the nuclear 
interaction only. The Coulomb interaction between nuclei and the X° particle is not included 
since we assume that the X° has a zero charge. The potential is supposed to be spherically 
symmetric. We solve the two-body Shrodinger equation by a variational calculation using 



the Gaussian expansion method 90|, and obtain binding energies. 



The two-body Shrodinger equation for a spherically-symmetric system is 

("C"^' + Vir) - i,{T) = , (11) 

where h is Planck's constant, fi is the reduced mass, V{r) is the central potential at r, E is 
the energy, and tp{r) is the wave function at r. If the mass of the X^ particle, i.e., mx, is 
much heavier than the light nuclides, fi is approximately given by the mass of the nuclide. 

The adopted RMS nuclear matter radii and their references are listed in columns 2 and 
3 in Table [H Binding energies of ground state X-nuclei are calculated with the interaction 
strength 6 (6^) and the mass mx taken as parameters. The obtained binding energies are 
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X Mass mx (GeV) 

FIG. 1: Contours of binding energies between nuclei and an corresponding to 0.1 MeV for the 
case of the Gaussian XN potential. Numbers attached to the contours indicate mass numbers of 
nuclei. 

used for the estimation of Q-values of various reactions as described below. Similarly, we 
calculate binding energies of nuclear excited states of ^B.e*x and ^Be^c with relative angular 
momentum of L = 1 by solving equation ( !TT|) for the L = 1 states. 

Figure [1] shows the contours of binding energies of ground state X-nuclei from mass 
number A = 1 to 9 in the case of the Gaussian type XN potential. The contours correspond 
to binding energy (BE) of BE = 0.1 MeV. This value of energy is chosen since weakly bound 
X-nuclei of BE < 0(0.1 MeV) tend to be photodisintegrated by background radiations 
during BBN epoch. Such weakly bound nuclei can not attain their large abundances without 
suffering from destruction processes. In a parameter region located at right upper side from 
a contour, the X-nucleus can form during BBN epoch, and thus possibly affects BBN. 

Figure [2] shows similar contours of binding energies of BE= 0.1 MeV in the case of the 
well type XN potential. The shapes of contours in both potential cases are very similar and 
change slightly. 

We adopt the Gaussian XN potential in calculating reaction rates and performing a 
network calculation of X-catalyzed BBN. After we introduce a mechanism of ^Be destruc- 
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TABLE I: Binding energies of X particles to Nuclei 
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^Derived by {r^^^f = {r^^^f - {a™^^f with a™^ = 0.875 ± 0.007 fm using a RMS proton matter 

radius determined in experiment as a RMS charge radius. 
"^Taken from ^He radius. 
''Taken from ^Li radius. 
'^Taken from ^Be radius. 
^Taken from radius. 



tion (this section) and show a result of the nucleosynthesis for the Gaussian XN potential 
(Sees. IIII AI and IIIIB|) . we delineate parameter regions for the ^Be destruction of not only 
the Gaussian but also the well XN potentials in Sec. IIII CI 
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FIG. 2: Contours of binding energies between nuclei and an corresponding to 0.1 MeV for the 
case of the square well XN potential. Numbers attached to the contours indicate mass numbers 
of nuclei. 



B. Reaction Rates 



Thermonuclear reaction rates {av) are roughly written (e.g., [95|, |96|) as 

(^^')nR = j2lklT)y^ r Ea{E) eM-E/ksT) dE, (12) 

where a is the cross section, v is the relative velocity, fi is the reduced mass of the system, 
ks is the Boltzmann constant, T is the temperature, and E is the kinetic energy in the 
center of mass system. Since there is a relation, i.e., E = the equation is identical to 

(^^)nr = ^^^^^3/, I a{E)vexp{-E/kBT)E'/' dE. (13) 

Reactions of a neutral particle and charged nuclei occur without the effect of Coulomb 
repulsion. If the quantity, i.e., a{E)v, is approximately given hj a + bE a.s a. linear function 
of E, then the reaction rate is simply given by 

Na{ctv)^r = NA(a + hksT) (14) 



2 



10 



where the Avogadro's number A^^ = 6.022 x 10^^ was multiphed to both sides of the equation. 
When the product, a{E)v, does not change drastically around some fixed point of E, the 
integral in equation (1131) receives a contribution from an energy region below E ~ ksT. 
Information from cross sections in higher energies is, therefore, not involved in the integral. 

Reactions triggered by two charged particles are, on the other hand affected by the 
Coulomb force. The astrophysical S'-factor is defined as 

5(E)=a(E)Eexp(Ml^^, (15) 

where Zi and Z2 are the charge numbers of interacting particles. The S'-factor might be well 
described by a linear function of E, i.e., S{E) = S{Eq) + a{E — Eq) = S{0) + aE, where 
Eq = 0.122{ZiZ2y^^A^^^Tg^^ MeV is the most effective energy in the integral in equation 
( |T2l) . We defined A = fi/{l amu) and Tg = T/(10^ K). The reaction rate can then be 



written 15|] in the form of 



A^aMnr = 7.82 X lOM Tg'/'exp ^12; 



/i / V ke Vb 



( (a/h) 



l22{ZlZlAf/^Tl'^ + ri.STgl i cm^ s-^moP^ie) 



[S(0)/(keVb)] 

where the Avogadro's number A^^ was multiplied, and lb=10~2^ cm^ was used. 

We estimate rates of several important reactions in this study. For both of reactions by 
a neutral plus a charged particles and those by two charged particles, calculated rates are 
used to derive linear fitting functions to be adopted at the energy range relevant to BBN, 
i.e., Tg < 1. 

We here assume that the mass of the X° particle, i.e., nix is 100 GeV. Calculations of 
nucleosynthesis are performed assuming the Gaussian XN nuclear potential as set up in 
Sec. Ill A[ We show results of the BBN catalyzed by the particle for two cases of different 
strengths of XA^ interaction, i.e., 5 = 0.1 (Case 1) and 0.2 (Case 2) in what follows. 

In this scenario of BBN catalyzed by the X^ particle, the ''Be can be destroyed at its X- 
capture (Sec. IllB lb ). The efficiency of this destruction, however, depends upon the fraction 
of the X^ particle escaping from the capture by ^He (Sec. IIIB 2b ). Since other reactions of 
X-nuclei can lead to productions of heavy nuclei, rates for such reactions are also estimated. 
Using reaction rates estimated as described in Sees. Ill BI L-4 in our nuclear reaction network 
(Sec. Ill B 51) . we perform a calculation of the catalyzed BBN. 
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The adopted reaction rates Na{(tv), per second per (mole cm~^), are shown in Tables HT] 
and mil Reaction Q- values are derived taking account of the calculated binding energies 
of the X-nuclei for Cases 1 and 2 (columns 4 and 5 in Table [T]). We use the notation, i.e., 
1(2,3)4 for a reaction 1 + 2 — )■ 3+4. Reaction rates related with the particle are estimated 
as follows. 



1. Non radiative reactions 



la. X(^Be,^He)^Hex 



This reaction is most important in this scenario. Its reaction rate is rather large since no 
radiation is involved in the reaction. It operates through the non resonant process unless 
there are resonant states lying near the energy level of initial scattering state. We then adopt 



the non resonant component of rate 97| for the normal nuclear reaction, i.e., ^Li(?7,, a)'^II as 
a rough approximation. The non resonant component of the rate for this reaction can be 
extracted most easily of all (n, a) reactions on light nuclides. This is because heavy nuclides 
have large densities of resonance so that they tend to have many resonant components for 



nuclear reactions 



The dependence of cross section on the reduced mass /i (or A in 



atomic mass units), i.e., a oc A has been used to correct for de Broglie wavelengths. 



lb. Other X-capture reactions 

Similarly, r eacti on rates of X(^Li, (i)^Hex, -^(^Li, t)^Hex, X(^Be,p)^Lix and 



X(^Be,n)^Bex 13l| are also taken from that of ^Li(fi,a)^H. They are corrected for the 
reduced masses. 

In Case 1, the Q- values of X(^Be,p)^Lix and XC^Be, ?T,)^Bex are negative. The reactions 
are then neglected as well as their inverse reactions which would never become important 
in changing the abundance of ^Be due to small abundances of ^Lix and ^Be^ as shown in 
Section TlTKl 
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Ic. Destruction of^Lix 

The main proton burning reaction of ^Lix is ^Lij5s:(p,^He)^Hex in Case 2. In Case 1, on 
the other hand, the Q- value of a three-body breakup reaction ^Li^ (p,^HeX)^He is positive. 
Its cross section is then much larger than that of ^Lix(p,^He)^Hex due to a larger phase 
space in the final state. This situation is the same as that in a similar reaction catalyzed 



by a long-lived negatively charged particle X~ [15|. The reaction rates of ^Lix(p,'^HeX)^He 
(Case 1) and ^Lix(p,^He)^Hex (Case 2) were both taken from that of ^Li(p, a)'^He. In the 
expression for a nonresonant contribution to the thermonuclear reaction rate [equation (fT6|) ]. 
the reduced mass was corrected. The survival of ^Lix thus differs from that of ^Li due to 
only the change in reduced masses of initial states relative to ^Li(p, a)'^He. 

Id. Production of ^ Bex 

Since the reaction ^Hex(«, 7)^Bex is found to be responsible for an accumulation of 
®Bex both in Cases 1 and 2 (see Sec. IIIip . flows of nuclear abundances toward higher 
mass numbers should be calculated. The most important reactions in this regard are 
^Bex{d,p)^Bex and ^Bex{d,n)^Bx The reaction rates of ^Bex{d,p)^Bex is taken from 
that of '^Bex{d,paYB.ex corrected for the reduced mass ^. The Q- values of ^Bex('^, n^Bx 
are negative. This reaction can, therefore, be neglected in environments of relatively low 
temperatures such as BEN. 

2. Radiative nuclear reactions 

We estimate rates of radiative capture reactions which can be important to leave signa- 
tures of the X° particle on primordial abundances of light elements. Wave functions of bound 
and continuum states are derived with the code RADCAP published by Bertulani 



which 



was modified and given proper input parameters described below. The code also calculates 
cross sections of forward and reverse reactions, and astrophysical ^-factors. 
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2a. X (a, -if Hex 



This reaction is very important for the existence of parameter region in which the pri- 
mordial ^Li abundance is reduced. The abundance of the free X° decreases if the strongly 
interacting X° particle is quickly captured by '^He after ^He nuclei are produced abundantly 
in the BBN epoch. Heavy nuclei such as ^'^Li and ^Be whose abundances build up after the 
'^He production are then not affected by free X° particles. The ^He^ nuclei can not react 
efficiently with Li and Be nuclei due to large Coulomb repulsion forces by them. 

Note also that the reaction between ^He and is mainly through a radiative capture. 
The critical binding energies of and ^He to an X° at which Q-values of X{a,NYAx 
reactions change from negative to positive are very high (see Table HVTl . These high binding 
energies are not realized in the cases of relatively weak XN interaction considered in this 
investigation. 

Whether a given reaction operates efficiently enough to change an abundance of relevant 
particle species is roughly determined by a comparison of its rate F and the cosmic expansion 
rate H. For a reduction of the X^ abundance via a radiative capture by ^He, the T / H ratio 
is given by 

I = (X\ ( ^ ^ ( ^ \ I 

H V0.25/ V6.2 X 10-10/ VO.l MeV/ V1.9 x 10^ cm^ s-i mol" 

where Y is the ratio by mass of ^He to total baryon, and r] is the baryon-to-photon ratio. 
\iT / H > 1, then the reaction X(a,7)^Hex quickly decreases the abundance of the free X° 
particle. 

In Case 1 the ground state of ^He^ exists although excited states do not. Due to the 
selection rule for electromagnetic multipole transitions, an electric dipole (El) transition 
from a relative s-wave scattering state into the ground s-wave bound state is not allowed. 
The cross section then has a predominant contribution from an El transition from a p-wave 
scattering state. 

In Case 2 there is one excited state of ^Hex with relative angular momentum of L = 1 (see 
Table [T]). The selection rule then allows an El transition from a relative s-wave scattering 
state into the excited p-wave bound state. 

For both cases the nuclear potential is given by equation ( !T0|) in calculating wave functions 
of the ground and excited states and scattering states of the ^He and X^ system with the code 
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RADCAP. Adopting = 1.59 fm, i.e., 6 = 1.30 fm, the potential is given numerically by 

V^{r) = -123 MeV5 exp {-[r/(1.97 fm)]^} . (18) 

2b. ^Hex(d,-ifLix 

We calculate the radiative capture cross section with the two-body model for the system 
of a ^Hex and a deuteron. The cross section is proportional to the electromagnetic matrix 



element squared [99|]. The matrix element is an integral over space of the scattering and 
bound state wave functions and the operator. The El operator is estimated from the initial 
state wave function. The operator is eirY'i^(f) for the radius r from the center of mass 
with the effective charge ei = einnYiex ~ '^'"^d) / (jn'^ncx + "^d)' ^^"^ ei — )■ e in the limit 
of an infinitely massive particle. Although the wave function of the ^Hex^c? system 
should be used in calculating the matrix element, as one method we approximately take 
the wave function of the ®Li-X bound state calculated with equation f|T9l) as that of the 
^Yiex^d bound state. As described hereinbelow, we try another method in which the wave 
function of the ^Lix nucleus is calculated by the two-body model for the ^Hex^c? system 
with its interaction tuned to reproduce binding energies of ®Lix. Three-body (^He, d and 
X^) quantum mechanical calculations are necessary for more consistent estimations of cross 
sections without the approximation. 

The nuclear potential for the ground state of ^Lix in the final state is given by equation 
fITU]) in calculating the bound state wave function only. Adopting = 2.35 fm, i.e., 

b = 1.92 fm, the potential is given numerically by 

V^{r) = -99.2 MeV5exp {-[r/(2.43 fm)]^} . (19) 

The nuclear potential for initial scattering states of '^Hex and d was taken from that 
between X° and d and that between ^He which is bound to the X° and d. The former is 
given by equation ( ITOl) . Adopting = 1.97 fm, i.e., b = 1.61 fm, the potential is given 
numerically by 

y/-'^(r) = -45.0 MeV(5exp {-[r/(2.19 fm)]^} . (20) 

The latter is given by 



V{r)T^-' = I p{r')Vr\x) dr' 
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27r /"^ 



dr' r'p{r') 



r+r' 



dx a;K 



a—d I 
N I 



X) 



(21) 



r Jo J|r-r'| 

where r is the radius from the center of mass of ^Hex to that of rf, r' is the distance between 
the position of ^Hex to that of the "^He nucleus, x = r + r' is the distance between the 
deuteron and ^He, and p{r') is the distribution function of ^He, V^~'^{x) is the nuclear 
potential between a free a particle and a deuteron, and spherical symmetries in p{r') and 
V-^~'^{x) were assumed in the last equality. 

The potential V^"*^ is given by the two-component Gaussian function 



a—d/ 



(22) 



1=1 



where Vi = 500 MeV, n = 0.9 rm, V2 = -64.06 MeV, and rs = 2.0 fm [15[. The potential 
V{r) ^ is then rewritten in the form of 



.a{B)-d 

In 



dr' r' p{r') I exp 



exp 



• (23) 



The Coulomb potentials for initial scattering states of ^Hex and d originate only from 
that between ^He (inside "^Hex) and d. It is approximately given [ist by 



Zde^ 



(24) 



where Zi^^ex = 2 is the electric charge of ^Hex nucleus, 64He^ and hd are ranges for charges 
of '^Hex and deuteron, respectively, hd = 1.47 fm is assumed which is derived from bd = 



2/3rJ^^^ and the RMS charge radius value of rg^^ = 1.4696 fm (lOO|. We take as the 64 
value the radius at which the charge density pc('") of ^Hex is exp(— 2) times the maximum 
value, i.e, pc(0). 

The charge density of Ax is given by 



Pc,Ax {r) 



PAAr')PcAr")dr' 
2Ze f°° 

dr' r'pAx (^') { exp 



nbr "'0 



— exp 



1^ I 



(25) 



where r' is the radius from the center of mass of Ax to that of A, r" is the position vector of 
charge contributed by the nucleus, r = |r"' + r"| is the distance from the center of mass of Ax, 
pyi^(r') is the distribution function of A, and pc,A{f") is the charge density of the nucleus. 
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The nuclear charge density was assumed to be pc,A{f") = ^e(7r6^)~^/^ exp [— in the 
last equality. The range parameter, i.e., 6 for ^He is given by 64 He =^j2/3rQ^^ and the RMS 
charge radius value of r^^^ = 1.80 fm 93 1. 



The charge densities of ^Hex are calculated, and the ranges of charge distribution are 
derived: 14^^^ = 2.97 fm {6 = 0.1) and 2.64 fm {6 = 0.2). They are used in equation (24). 

Figure [3] shows the adopted nuclear potential V^~'^ and the calculated distribution func- 
tion, i.e., p, and charge density, i.e., pc, of ^He^, and nuclear {V^^^''~'^) and Coulomb 
(V^^^-*"*^) potentials (solid lines) for Case 1 (5 = 0.1) as a function of radius from the center 
of mass of ^Hex- The center of mass is approximated to be the position of since the 
mass of X° is much larger than that of nucleon. 

Figure m shows similar potentials and densities (solid lines) for Case 2 {6 = 0.2). 

We calculate rates using another method for the same reaction. We assume that the ^Lix 
nucleus in the final state is described as the two-body bound state of the ^Hex^c? system. 
Both of initial and final state wave functions are generated from a ^He-d potential tuned 
to reproduce the binding energy of ^Lix relative to separated ^Hex and deuteron consistent 
with results in Table [T] [2.44 MeV for Case 1 {6 = 0.1) and 2.16 MeV for Case 2 (5 = 0.2)]. 
We use potential terms given by equations (j20l) . ( !23l) and (j24l) . The parameter 6 in the 
nuclear potential between an X° and a deuteron [equation ( l20l) ] is fitted: 6tit = 0.494 for 
Case 1 and 6fit = 0.373 for Case 2. Calculated reaction rates are as follows: 

NA{av) = 1.9 X 10%"^/^ exp(-8.45/r9/^)(l + 2.6X9^^ + O.78T9) cm^ s'^moP^ (26) 

for Case 1 (5 = 0.1), and 

Na{(tv) = 1.7 X 10^T9"^/Sxp(-8.45/T9/^)(l + 3.6r9/^ + l.ir9) cm^ s'^moP^ (27) 

for Case 2 {6 = 0.2). These rates differ from our standard rates in Tables HTl and lllll by 17-26 
% (Case 1) and -45- -33 % (Case 2), respectively, in the important temperature range of 
T9 =0.5-1 corresponding to Eq =0.15-0.25 MeV. 

2c. '^Hex(a,'yfBex 

In both of Cases 1 and 2 there is one excited state of ^Bex with L = 1 below the energy 
level of the initial separation channel of ^Hex and a. The excitation energies are 3.72 MeV 
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2 4 6 

r (fm) 



FIG. 3: Nuclear potentials between an a and a deuteron {V^ '^), an and a deuteron (y^~'^), 
and the Coulomb potential between an a and a deuteron (Vq as a function of distance from 

the center of mass of the ^Hex+c? system (solid lines). The similar potentials for the ^Hex+o: 
system (dashed lines), the distribution function (p), and the charge density (pc) of ^Hex (solid 
lines) are also drawn. It was assumed that the mass of particle is mx = 100 GeV, and that 
the interaction strength of XN force is 0.1 times that of NN force {5 = 0.1). 



{6 = 0.1) and 5.89 MeV {6 = 0.2), respectively. Both rates for final states of the ®Bex 
ground and excited states are then calculated. 

The nuclear potential for the ground state of ®Bex in the final state is given by equation 
f IlOp in calculating the wave function. Adopting = 2.33 fm, i.e., b = 1.90 fm, the 

potential is given numerically by 

V^{r) = -134 MeVSexp {-[r/(2.41 fm)]^} . (28) 

The nuclear potential for initial scattering states of ^He y and a was approximately taken 
from that between two a particles given by the three-component Gaussian function 

3 

V^{r) = Y^Viexp[-{r/r,f], (29) 

i=l 
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FIG. 4: Same as in Figure [3] when the interaction strength of XN force is 0.2 times that of NN 
force {b = 0.2). 



where Vx = -1.742 MeV, n = 3.00 rm, V2 = -395.9 MeV, = 1.90 rm, V3 = 299.4 MeV, 



isl- 



and r3 = 1.74 fm 

The adopted nuclear potential V^~°' and the calculated nuclear (V^^^^ ") and Coulomb 
{Vq^^^ ") potentials (dashed lines) for Case 1 (5 = 0.1) are shown as a function of radius 
from the center of mass of ^He^ in Figure El Similar potentials (dashed lines) for Case 2 
{6 = 0.2) are shown in Figure HI 

We calculate cross sections by the two-body model in the same procedure as performed 
in Sec. Ill Bb 2. First and second terms of reaction rates in Tables ITTl and UTTl are calculated 
rates for transitions to the ground and excited states, respectively. The reaction leading 
to the ground state is via an El transition from an initial p-wave scattering state, while 
that leading to the excited state is predominantly an El transition from an initial s-wave 
scattering state. 
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TABLE II: Reaction rates for {mx, 6)={100 GeV, 0.1) 



Reaction 


Reaction Rate (cm^ s ^ mole ^) 




Energy (MeV)" 


Reverse Coefficient'' 




X (6 Li, d)-* Hex 


3.8 X 10" 






4.919 


14.556 


'^Ux{p,^Hea)X 


3.6 X 10iOTg"^''^exp(-8.81/Tg^^) 




— 




35.421 




1.9 








44.225 




23xin^-l-2fix1 n^To 






7.472 


31.659 


4Hex(rf,7)'^Lix 


1.0 x lO'^Tg"^''^ cxp(-8.45/Tg/^)(l + 6.9T^^^ + 


2.OT9) 


0.15-0.25 


2.717 


28.331 


X(7Li,t)'iHex 


2.9 X 10*5 






6.748 


3.034 


X(''Be,3He)4Hex 


2.9 X 10" 






6.748 


13.252 


■*Hejf (q,7)»Box 


7.8 X 105Tg"^/^exp{-16.79/Tg''^)(l - 0.39Tg/^ - 
+4.8 X lO^Tg"^/^ exp(-16.79/Tg''^)(-l + 2.1Tg''^ 


O.O58T9) 
+ O.32T9) 


0-1 
0.3-0.5 


7.487 


45.509 


8Bex(rf,p)«Bex 


9.8 X 10"Tg"^/^exp(-13.49/Tg''^)'' 






1.060 


7.242 



"Rates are estimated by linear fits of S'-factors in the energy ranges. 

''For nuclides a = i,j,k,... with mass numbers Aa and numbers of magnetic substates ga, the reverse 



coefRcients are defined as in Ref. 



10^: 0.9867{g,gj/gk){A^A,/Ak^f/^ for the process ix{j,'j)kx, and 



{9i9j/i9kgi)][Aij,Aj/{AkAi^)f/^ for the process ix{j,k)lx- 
Qg = 11.605(Q/M|;V). 

J] is taken. 



The rate in Ref. 



3. P- decay 

The reaction rate of ^Bex{,P~^'^e)^^ix is estimated using that of ^He(,/3~z7g)^Li. It was 
corrected for a phase space factor related to the reaction Q- value. The rate is then given by 
Tx = (ln2/Ti/2)(Qx/<5)^, where T1/2 is the half life of ^He, Qx and Q are the Q- values for 
the /3-de cay o f ^ Bey and ^He, respectively. T1/2 = (806.7 ± 1.5) ms and Q = 3.508 MeV are 



adopted 



101| . See Tables [TTl and [TTIl for derived Qx-values. 



4- Unimportant pathways through ^Hex and ^Lix 

Reactions of ^Hex and ^Lix are found to be unimportant in this study under the as- 
sumption of our models described above. 

The binding of the particle to ^He and ^Li can lead to stabilizations of such bound 



states against neutron and proton emissions, respectively ^ | . This study, however, indicates 
that ^Ax nuclei do not play a significant role in the BBN epoch in both Cases 1 and 2. The 
reason is as follows: 

A sufficient condition for that the ^Ax nuclei are not important for a production of 
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TABLE III: Reaction rates for {mx, 5)=(100 GeV, 0.2) 



Reaction 


Reaction Rate (cm^ s ^ mole ^) 




Energy (MeV)" 


Reverse Coefficient'' 




X (6 Li, d)-* Hex 


3.8 X 10'' 






4.919 


96.708 


XCBe,nfBex 


2.9 X 10'' 




— 


34.139 


0.130 


XCBe,pY'Ux 


2.9 X 10" 






11.380 


56.727 


6Lix(p,3He)4Hex 


3.6 X 10iOTg"^''^cxp(-8.81/r9''^) 




— 


0.593 


38.678 


6Bex(,e+i/e)«Lix 


6.2 X IQ-io 








44.225 




■SOxlO'' — 22x lO'^Tn 






7.472 


113.811 


''Hex(rf,7)''Lix 


8.2 x 103Tg"^''^exp(-8.45/Tg''^)(-H-6.6rg''^ - 


1- 2.OT9) 


0.15-0.25 


2.717 


25.074 


XC^Li,t)4Hex 


2.9 x ID" 






6.748 


85.186 


X(''Be,3He)-4Hex 


2.9 X 10'' 






6.748 


95.404 


4Hex(«,7)*Bex 


2.4 X lO'^Tg'^''^ exp{-16.79/Tg''^)(l + CIST^^^ + 
+2.7 X lO^'Tg^^''^ exp(-16.79/T9''^)(-l + 2.7Tg''^ 


O.O22T9) 
+ O.4IT9) 


0-1 
0.3-0.5 


7.487 


82.366 


8Bex(rf,p)SBex 


9.8 X 10"Tg"^''^exp(-13.49/Tg''^)'' 






1.060 


20.957 



"Rates are estimated by linear fits of S'-factors in the energy ranges. 

''For nuclides a = i,j,k,... with mass numbers Aa and numbers of magnetic substates ga, the reverse 



coeflticients are defined as in Ref. 



10^: 0.9867(5,ffj/5fc)(^,_^^,Mfc^j3/2 for the process ix{j,'y)kx, and 



[9i9j / i9k9i)][Ax / i^kAx)?^'^ for the process ix{j,k)lx- 
'Qg = 11.605(Q/Mi!V). 



^'The rate in Ref. 



J] is taken. 



heavy nuclei of A > 6 is satisfied in Cases 1 and 2: Q-values of reactions ^Hex('i, nYLix and 
'^B.ex{d,p)^B.ex are negative. This means that ^Li^ and ^He^ produced during BBN (if any 
via radiative captures of p or n by ^Hex) is quickly processed by n or p back into ^Hex via 
^Lix{n,d)^B.ex and ^B.ex{p,d)^B.ex- Note that the destruction reactions are strong since 
no radiation is involved in the reactions. 

In passing, if the energy level of ^Hex (^Lix) nucleus is lower than that of ^Hex+^ (p) 
separation channel, i.e., Q(*B.ex + N — J-^ Ax) ^ MeV, the X-nucleus can be produced 
during BBN. Amounts of such X- nuclei are, however, very small due to the reason explained 
in the previous paragraph. In this investigation the ^Ax nuclides have not been stabilized 
in both cases of potential types for parameter space studied. The ^Ax nuclei then do not 
play a role. There is, however, a large uncertainty in binding energies of X-nuclei which 
stems from XA potentials, and from precise nucleon density of X-nuclei. Such uncertain 
points should be studied using specific particle models describing potentials and dedicated 
quantum many-body models taking account of interactions inside the normal nuclei such as 
^A =^He+A^ during the processes. 
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5. Nuclear Reaction Network 



We calculate binding energies of X-nuclei, and Q-values of reactions involving the X° 
particle and X-nuclei. We estimate rates of several (forward and reverse) reactions which 
possibly play roles in BBN as described above, and calculated the catalyzed BBN. We note 
that productions of double X-nuclei, i.e., Axx are not taken into account in this study. The 
ratio of the rate for reaction between X-nuclei and an to that of cosmic expansion is 
given by 



r 
H 



'nx/ub 
10-4 



T 



Na{(tv) 



6.2 X 10-10/ VO.l MeV/ \ 1.2 x 10^ cm3 s"! mol" 



(30) 



An efficient operation of a reaction should satisfy T/H > 1, which is possible when Nj^{av) ^ 
10® cm^ s-i mol-i. X-nuclei with significant abundances during the catalyzed BBN seen in 
Sec. mil are ^Hex, ^Lix, ^Bex and ^Bex- ^Hex can react with an X" particle only via the 
radiative capture since there is no exit channel of particle break up. Since reaction rates 
of radiative capture are typically small, ^B.exx production is not efficient. The effect of 
^Hexx would thus be negligible although some fraction of ^Hex would react with an X^ and 
form the ^Hexx- Other X-nuclei, i.e., ®Lix, ^Bex and ^Bex can, on the other hand, react 
nonradiatively with an X^. Possible reactions are ^Lix{X,dYB.exx, '^Bex(X, a)4Hexx, 
®Bex(X, ax)^Hex, ^Bex(X, n)®Bexx and ^Bex(X, Q;n)4Hexx (whether their Q-values are 
positive or negative should be determined by more sophisticated estimations of binding 
energies). Since abundances of ®Lix, ^Bex and ^Bex are rather small in cases studied in 
this paper (Sec. IIIip . their processing would have negligible effects on final abundances of 
light elements. We note, however, that final abundances of ®Lix and ^Bex may be reduced 
through nonradiative reactions with an X" pa rticle. 



The BBN network code of Refs. 



103 



104| is modified and used. The X" particles and 



relevant X-nuclei are included as new species. Reactions connecting normal and X-nuclei 
and the X^ particle are added to the co de (see Tables HTl and IITTI for their rates). Nuclear 



reaction rates for the SBBN 
Refs. 58 



103 



Il04j have been replaced w ith n ew rates published in 



105[ | and the adopted neutron lifetime is r„ = 881.9 s [106 1. 
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FIG. 5: Calculated abundances of normal and X-nuclei (solid lines) as a function of Tg. The 
mass of particle was assumed to be 7nx = 100 GeV, and the interaction strength of XN 
force is 0.1 times that of A^A^ force {5 = 0.1). For this figure, we took the X^ abundance to be 
Nx/nb = 1.7 X 10~^ {Yx = Nx/s = 1.5 x 10"-^^), and its lifetime to be much longer than BBN 
time scale, i.e., tx ^ 200 s. The X^ reaction rates are given as described in the text. Section HIl 
The dashed lines correspond to the abundances of normal nuclei in the standard BBN model. 



III. RESULTS 



A. Evolution of nuclear abundances 



Figure |5] shows results of the abundances of normal and X-nuclei as a function of tem- 
perature for Case 1 (5 = 0.1). The mass of the X° has been taken to be mx = 100 GeV. 
Its initial abundance is Nx/nt = 1.7 x 10~^ (Yx = Nx/s = 1.5 x 10~^^), where Nx and rib 
are the number densities of the particles and baryons, respectively. This abundance is 
chosen as an example leading to a significant ^Be reduction. It is set as a parameter here 
since the abundance is very uncertain. The decay lifetime is assumed to be much longer 
than BBN time scale, i.e., tx ^ 200 s so that effects of the decay are not seen. The X° 
particle is assumed to have been long extinct by now. 
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FIG. 6: Same as in Figure 3 when the interaction strength of XN force is 0.2 times that of NN 
force {5 = 0.2). 



At high temperatures Tg > 1, the particles exist mainly in the free state since efficient 
photodisintegrations of X-nuclei destroy the bound state. At Tg ~ 1 the ^He synthesis occurs 
as in SBBN, and about one third of X° particles are then captured by ^He nuclei. ^Hex 
nuclei produce d in this epoch react with normal nuclei, and affect abundances of ^Li [by 
^Hex(tjLi)X I132I], ^Lix [by ^Hex(ci, 7)^Lix], and ^Be.Y [by ^Hex(a, 7)^Bex]. Note that 
^Lix nuclei produced at Tg ~ 1 experience a strong destruction process, i.e., ^Li(p,^Hea)X. 
^Bex is produced by ^Bex{d,p)^Bex- At last, the most important processes operate. Free 
particles which survived the capture by ^He react with ^Be [by Ar(^Be,^He)^Hej!s:] and 
''Li [by X(^Li, t)^Hej5s:]- The abundances of '^Be and ^Li thus decrease. 

Figure |6] shows results of the abundances of normal and X-nuclei as a function of tem- 
perature for Case 2 {6 = 0.2). Parameters other than S are the same as in Figure [5l A clear 
difference from Case 1 is a complete capture of the X^ particle by ^He (see Section [11 B 2^ . 
Decreases in the abundances of ^Be and ^Li are, therefore, not seen. The productions of 
^Lix, ^Bex, and ^Bex result following the ^Hex production similarly to Case 1. 
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B. Decrease in the primordial Li abundance 



Figure [7] shows abundances (solid curves) of ^He (mass fraction), D, ^He, ^Li and ^Li 
(by number relative to H) as a function of the baryon-to-photon ratio t] or the baryon 
energy density ^2^^^ of the universe. The solid curves are the calculated result in the X 
catalyzed BBN for Case 1, i.e., (m^, S, Yx, rx) = (100 GeV, 0.1, 1.5 x 10-^^ oo). The 
dashed curves are those in the SBBN. The boxes correspond to the adopted constraints 
on primordial abundances (see Appendix ^. The vertical stripe shows the 2a limits on 



fl^h^ = 0.02258;{|;S provided by WMAP |59| for the ACDM+SZ+lens model. 

The decrease in the ''Li abundance is found, while other nuclear abundances are not 
changed. A solution to the ^Li problem is thus found in this model. One should note, 
however, that the abundance of ^Bex could be higher than the adopted constraint on the 
primordial ^Be abundance. For example, in Figure |5l the final abundance of ^Be^/H ~ 
2 X 10~^^ is shown. Effects of the decay of the particle inside X-nuclei are not addressed 
in this paper. They should be studied in order to estimate fractions of ^Be^ which remain 
as ^Be after the decay of X^ particle. 

C. Parameter region for Li reduction 

If the strength of XN interaction is relatively weak as in the Cases 1 and 2 which we 
study in this paper, most strong reactions for the X^ particle to get bound to nuclei would 
be non-radiative X°-capture reactions which are found important in the present model. We, 
however, note that efficiencies of the reactions are uncertain, and that the present result is 
based on a rough assumption that the non-radiative cross sections have been set equal to 
the ^Li(n, a)^B. cross section excepting for the factor of reduced mass. Although effects of 
heavy X^ particles during BBN was first studied recently ^ L the possibility of X^ capture 
reactions via nucleon emission has been mentioned in 1995 2[- In order for these reactions 
to occur efficiently during BBN, reaction Q-values need to be positive. The Q-value of the 
reaction X{A, b)Bx is given by 

g = BE(fi) + BE(5x) + BE(6) -BE(A), (31) 

where BE(A), BE{B) and BE(6) are the binding energies of nuclei A, B and h with respect 
to separated nucleons, respectively, and BE(i?x) is the binding energies of Bx with respect 
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FIG. 7: Abundances of ^He (mass fraction), D, ^He, ^Li and ^Li (by number relative to H) as 
a function of the baryon-to-photon ratio t] or the baryon energy density Qsh^ of the universe. 
The sohd curves are the calculated results in the X catalyzed BBN for the case of {mx, (5, Yx, 
rx)=(100 GeV, 0.1, 1.5 x 10~^^, oo), while the dashed curves are those in the standard BBN. There 
is virtually no difference between the dashed and solid cur ves f or ^H e, D, ^He , and ^Li. T he b oxes 
represent the adopted jibundance constraints from Refs. 
■^He, 5J] for ^Li, and [52] for ^Li, respectively. The vertical stripe represents the 2 a ri^/i^ limits 



1071 . 



IMi for ^He, 



1091 ] for D, [llfi] for 



provided by WMAP [5£ 
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TABLE IV: Critical binding energies of X-nuclei 



IVniTil^pr 


part ion 


RinHinp" ptipvp'v (TVTpVi 


1 

X 


v\ n, 11) XIX 


9 994 




£i 




9 994 


\j 


w 11, iix 


fi 957 


A 




5 4Q4 


5 


lie, jJ) iix 


19.815 






90 578 


7 




4 4^7 


Q 

o 




5 "^Q 


q 




Q Q75 


If) 




7 950 


1 1 

± ± 




5 fiDfi 

fj . \J\J\J 


12 
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to separated B and X". 

In Table HVl critical binding energies of Bx realizing Q > are listed. We have tak en d ata 
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on binding energies of normal nuclei from TUNL Nuclear Data Evaluation Project 

More important reactions than nucleon emissions in th is X catalyzed BBN scenario are 
X(6Li,d)^Hex, X(^Li,t)^Hex, and X(^Be,3He)^Hex l3- Critical binding energies realiz- 
ing Q > for emissions of particles other than nucleons are also listed in Table llVi 

Figure E] shows contours in the parameter space (mx, 5) for critical binding energies of 
X-nuclei (thin and thick smooth curves) for the case of the Gaussian XN potential. Critical 
binding energies chosen for the plot are those of reactions, whose numbers are defined in 
TablelVl 2 (for n^), 1 (^H^) 3, 4 (^Hx), 5 (^Hx), 6 (^Hex), 14, 15, 16 (^Hex), 7 (^Hex), 8 
(^Lix), 10, 11 (^Lix), 12 (^Bex)- Numbers attached to the contours indicate mass numbers 
of nuclides for elements which have more than two isotopes plotted. The contour of Q = 
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0.1 1 10 100 



X Mass (GeV) 

FIG. 8: Contours of binding energies between nuclei and an particle corresponding to Q = 
of reactions (thin and thick smooth curves) for the case of the Gaussian XN potential. Numbers 
attached to the contours indicate mass numbers of nuclei. Above the contours, reaction Q-values 
are positive (see text, Sec IIII CJj) . Zigzag curves correspond to boundaries above which an excited 
state of ^He* (upper line) and ^Be* (lower) exist, respectively. In the shaded region, the ^Li problem 
can be resolved. 



for the proton decay of ^Bex, i.e., ^Bex{, 2p)^Hex, is also shown as a thick solid line. Above 
the contours, reaction Q-values are positive. Zigzag curves correspond to boundaries above 
which an excited state of ^He* with L = 1 (upper line) and ®Be* with L = 1 (lower) exist, 
respectively. 

Figure [9] shows contours for the case of the square well XN potential corresponding to the 
same boundaries as in Figured The contours are very similar to those in Figure |8] excepting 
for that of the proton decay of ^Bex- Although the difference in the contours of ^Bex does 
not affect the BBN significantly, a realistic estimation of binding energies regarding all light 
nuclei are very important to obtain a realistic result of the catalyzed BBN. 

We delineate the parameter region which might be responsible for a reduction of the 
primordial ^Li abundance. The baryon-to-photon ratio is fixed to he rj = 6.225 x 10~^° from 



WMAP determination 



59| . In both Figures [8] and [9l it is found that the contours of the 
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FIG. 9: Same as in Figure [8] for the case of the square well XN potential. 



boundaries for existences of ^He3c are above the contours of the reaction X(^Be,^He)^Hex 
(the second lowest thin dashed lines). In the parameter region in right upper side from the 
curve of ^He^c, free particles are captured onto ^He at Tg ~ 1 before they can react with 
^Be to reduce its abundance. 

In the shaded region below that curve and above the curve of X(^Be,^He)^Hex, some 
amount of free X°s possibly remain, and they can reduce the ^Be abundance. This shaded 
region is, therefore, a possible parameter region where the ''Li problem is solved. For a 
significant destruction of ^Be, however, a relatively large value of initial abundance, i.e., 
Yx ~ O(10~^^) is needed. Under the assumption of thermal freeze-out from equilibrium of 
the particle, the initial (relic) abundance is large if the annihilation cross section of X° 
particle is small [c.f. equation The required abundance may then realize in the scenario 
of sub-strongly interacting particle X^ which has survived the annihilation due to its small 
interaction strength. 
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IV. CONCLUSIONS AND DISCUSSION 



We have investigated effects on BBN of a long-lived strongly interacting massive particle 
(SIMP) for different masses mx and strengths of XN interaction, i.e., 6. Binding energies 
of bound states of nuclei and an particle, i.e., X-nuclei, are calculated for two types of 
XN potentials, i.e., Gaussian and well types. It is shown that calculated binding energies 
are not largely dependent upon the potential shapes, and are determined by the interaction 
strength at a given mass of X^. 

Evolutions of light element abundances are calculated as a function of the temperature 
for two specific cases of relatively weak interaction strengths. We found that ''Be and ''Li can 
be destroyed by the nuclear capture reactions of free X° particles. The reactions identified 
as destruction processes are X('^Be,^He)^Hex and X(''Li, t)^Hex- We show that the lack 
of an excited state of ^Hex with a relative angular momentum L = 1 is essential for some 
fraction of the X^ particles to escape capture by ^He. 

We suggest that the ^Li problem could be solved based upon a net work calculation of 
catalyzed BBN, and found the parameter region in the {mx, S) plane where the ^Li problem 
can be fixed. 

We note that the results have been derived under the assumption that the X° particles 
do not change nuclear structure apart from sticking unaltered nuclei to the particles. This 
rough approximation is unlikely to be true especially in the case of relatively large strength 
of interaction since the XN potential is not much weaker than the NN potential and can 
not be neglected. More realistic estimations of wave functions and binding energies of X- 
nuclei need to include such changes in nuclear structures with the use of three or more-body 
models. In the case of strongly bound ^He, the binding energy of the ground state is 28 MeV 
which is larger than binding energy of ^He^ with respect to separated ^He and X^. The 
effect of change in nuclear density caused by the existence of the X^ then tends to be small 
for ^He. The calculated binding energy of ^Hex and the cross section for radiative capture of 
X° on ^He are, therefore, likely to be accurate. If the XN and XA potentials adopted in this 
paper describe well the real interaction, the main uncertainty in the BBN calculation would 
be in the estimations of nuclear reactions involving the X^ particle. In this paper, the cross 
section of the most important reaction, i.e., X('^Be,^He)^Hex was estimated only by analogy 
with ^Li(n, a)^H, and many other cross sections were also estimated using standard nuclear 
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reaction rates. Radiative reaction cross sections were estimated approximately within the 
framework of two-body models. The rates should, however, be calculated in more rigorous 
quantum many body models, not by the rough Born approximation, in order to derive 
realistic values. Although there might be errors in the calculated abundances in the X- 
catalyzed BBN of one order of magnitude or so, we argue at this moment that there is a 
possibility of ''Li reduction in the BBN model including a long-lived sub-SIMP. 

Effects of possible direct interactions of decay products of with the remaining nuclei 
A at the decay of X^ in an X- nucleus Ax are not taken into account yet. They should 
be studied in the future in order to better estimate final abundances of the light elements. 
In addition, nonthermal nucleosynthesis triggered by the decay process of the X^ particle 
might change the abundances of normal nuclei if the energy injection by the decay were large 
enough. Recent studies suggest that the radiative decay could lead to the production of ^Li 
to the level at most ~ 10 times larger than that observed in MPHSs when the decay life is 
of the order of ~ 10^ — 10^^ s which is associated with ^He production The hadronic 
decay, on the other hand, can be a solution of both the lithium problems although that case 



gives somewhat elevated deuterium abundances 35 



43|. 



For example, we assume that the mass and the initial abundance of the X are irtx = 
100 GeV and Yx = 1-5 x 10~^^, respectively. The energy injection at the decay of the X 
into hadronic jets is constrained to be ^ 0(1-100 GeV) if the lifetime is tx ^ 10^ s thr oug h 
abundances of D, ^He, ^Li and ^He depending upon the lifetime (figure 38 of Ref. 32[). 
The energy injection into electromagnetic particles is, on the other hand, constrained to be 
^0(10 GeV) if the lifetime is tx ^ 10^ s (figure 1 of Ref. j^^]). This amount of energy 
injection tends to attain ^Li production up to the observed level in MPHSs. 

We summarize a present status of several models which have effects and thus leave ob- 
servational signatures on primordial light element abundances. In the BBN model catalyzed 
by a long-lived sub-SIMP studied in this paper, the abundance of ''Li can be reduced below 
the level of SBBN prediction. In the BBN model catalyzed by a long-lived SIMP, the abun- 
dances of ^Be or B can be high j4|. Moreover, the isotopic abundance ratio, i.e., ^'^B/^^B can 
be high, which is never predicted in other scenarios for boron production {4]. In the BBN 
model catalyzed by a negatively charged massive particle (CHAMP), the ^Li abundance can 

lan (0.04 — 1) times as large as 



be high jsl. Only if the abundance of the CHAMP is more t 
that of baryon 12|, llll| . the ^Li reduction can be possible 



9[. A signature of CHAMP on 
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Be abundance has been estim ated to be negligible 12|, in the light of a rigorous quan- 



tum mechanical investigation 



by supernova explosion in an ear 



86 



112|. The cosmological cosmic ray nucleosynthesis triggered 
y epoch of the structure formation can produce ^Li 113 | 



as well as ^Be and ^"^'^^B 

of ^Be can be higher than in the SBBN |115l4ll9| . 



114j . In b aryon-inh omogeneous BBN models, the abundance 



Appendix: Constraints on primordial light element abundances 



We adopt constraints on primordial abundances as follows: 

A mean value for the primordial deuterium abundance in Lyman-a absorption syste ms in 
the foreground of high redshift QSOs has been estimated to be log(D/H)=— 4.55±0.03 109|. 
We adopt this value and a 2a uncertainty, i.e., 

2.45 X 10^^ < D/H < 3.24 x 10"^ (A.l) 

^He abundances are measured in Galactic HI I reg ions through the 8.665 GHz hyperfine 
transition of ^He^, i.e., '^He/H=(1.9±0.6) x 10~^ 110| . However, abundances in extragalactic 



objects have not been measured, and it is not known whet 
through the course of stellar and galactic evolution 120 



ler ^He has increased or decreased 



12 1| . ^He is more resistant to the 



stellar burning than deuterium. Because the deuterium abundance does not appear to have 
decreased since the BBN epoch until the solar system formation 122| . we do not assume 
a decrease in the ^He abundance after BBN in amounts larger than the uncertainty in the 
abundance determination of Galactic HII regions. Although a constraint on the primordial 
^He abundance is rather weak considering its uncertainty, we take a 2cr upper limit from 
abundances in Galactic HII region as a rough guide, i.e.. 



^He/H < 3.1 X 10" 



(A.2) 



For the p rimordial helium abunda nce w e adopt two different constraints, i.e, Y = 0.2565± 
0.0051 107| and y = 0.2561±0.0108 108| both from observations of metal-poor extragalactic 
HII regions. We take 2a limits of 



and 



0.2463 <Y < 0.2667 (ITIO), 



0.2345 <Y < 0.2777 (AOSIO). 



(A.3) 



(A.4) 
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An upper limit on the ^Li abundance is taken from the possible plateau abundance of 



6Li/H= (7.1±0.7) X 10~^2 observed in metal-poor halo stars (MPHSs) 
is included and we derive 

6t ; /xj ^ o K XX 1 n— 12 



52| . A 2 cr uncertainty 



^Li/H < 8.5 X 10" 



(A.5) 



A limit on the ^Li abundance is taken from observations of MPHSs, i.e., ^Li/H= 



{1.23+lil) X 10-1° (95% confidence limits) 
is then 



5J| . The adopted constraint on the ''Li abundance 



0.91 X 10"'" < 'Li/H < 1.91 X 10 



-10 



(A.6) 



Although we do not use constraints on abundances of nuclei with mass number A > 9 
in this study, a primordial ^Be may be related to the scenario (see Section IIIIBI) . An 
upper li mit o n ^Be abundance should be taken from the minimum abundances observed in 



MPHSs 



123|, 



I.e. 



^Be/H < 10 



-14 



(A.7) 



Acknowledgments 

This work is supported by Grant-in-Aid for JSPS Fellows No. 21. 681 7 (Kusakabe) and 
Grant-in- Aid for Scientific Research from the Ministry of Education, Science, Sports, and 
Culture (MEXT), Japan, No.22540267 and No.21111006 (Kawasaki) and also by World 
Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. 



[1] D. A. Dicus and V. L. Teplitz, Phys. Rev. Lett. 44, 218 (1980). 
[2] R. Plaga, Phys. Rev. D 51, 6504 (1995). 

[3] R. N. Mohapatra and V. L. Teplitz, Phys. Rev. Lett. 81, 3079 (1998), hep-ph/9804420. 
[4] M. Kusakabe, T. Kajino, T. Yoshida, and G. J. Mathews, Phys. Rev. D80, 103501 (2009), 
0906.3516. 

[5] M. Pospelov, Phys. Rev. Lett. 98, 231301 (2007), hep-ph/0605215. 
[6] K. Kohri and F. Takayama, Phys. Rev. D76, 063507 (2007), hep-ph/0605243. 
[7] R. H. Cyburt, J. R. Ellis, B. D. Fields, K. A. Olive, and V. C. Spanos, JCAP 0611, 014 (2006), 
astro-ph/0608562. 



33 



[8] K. Hamaguchi, T. Hatsuda, M. Kamimura, Y. Kino, and T. T. Yanagida, Phys. Lett. B650, 

268 (2007), hep-ph/0702274. 
[9] C. Bird, K. Koopmans, and M. Pospelov, Phys. Rev. D78, 083010 (2008), hep-ph/0703096. 
[10] M. Kusakabe, T. Kajino, R. N. Boyd, T. Yoshida, and G. J. Mathews, Phys. Rev. D76, 

121302 (2007), 0711.3854. 
[11] M. Kusakabe, T. Kajino, R. N. Boyd, T. Yoshida, and G. J. Mathews, Astrophys. J. 680, 846 

(2008), 0711.3858. 

[12] M. Kusakabe, T. Kajino, T. Yoshida, and G. J. Mathews, Phys. Rev. D81, 083521 (2010), 
1001.1410. 

[13] K. Jedamzik, Phys. Rev. D77, 063524 (2008), 0707.2070. 
[14] K. Jedamzik, JCAP 0803, 008 (2008), 0710.5153. 

[15] M. Kamimura, Y. Kino, and E. Hiyama, Prog. Theor. Phys. 121, 1059 (2009), 0809.4772. 
[16] M. Pospelov (2007), 0712.0647. 

[17] M. Kawasaki, K. Kohri, and T. Moroi, Phys. Lett. B649, 436 (2007), hep-ph/0703122. 
[18] T. Jittoh et al., Phys. Rev. D76, 125023 (2007), 0704.2914. 
[19] T. Jittoh et al., Phys. Rev. D78, 055007 (2008), 0805.3389. 
[20] T. Jittoh et al. (2010), 1001.1217. 

[21] M. Pospelov, J. Pradler, and F. D. Steffen, JCAP 0811, 020 (2008), 0807.4287. 
[22] M. Y. Khlopov and C. Kouvaris, Phys. Rev. D77, 065002 (2008), 0710.2189. 
[23] J. R. Ellis, D. V. Nanopoulos, and S. Sarkar, Nucl. Phys. B259, 175 (1985). 
[24] R. H. Cyburt, J. R. Ellis, B. D. Fields, and K. A. Olive, Phys. Rev. D67, 103521 (2003), 
astro-ph/0211258. 

[25] J. R. Ellis, K. A. Olive, and E. Vangioni, Phys. Lett. B619, 30 (2005), astro-ph/0503023. 
[26] N. Terasawa, M. Kawasaki, and K. Sato, Nucl. Phys. B302, 697 (1988). 
[27] M. Kawasaki et al., Nucl. Phys. B419, 105 (1994). 

[28] M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93, 879 (1995), hep-ph/9403364. 
[29] E. Holtmann, M. Kawasaki, and T. Moroi, Phys. Rev. Lett. 77, 3712 (1996), hep-ph/9603241. 
[30] M. Kawasaki, K. Kohri, and T. Moroi, Phys. Rev. D63, 103502 (2001), hep-ph/0012279. 
[31] M. Kawasaki, K. Kohri, and T. Moroi, Phys. Lett. B625, 7 (2005), astro-ph/0402490. 
[32] M. Kawasaki, K. Kohri, and T. Moroi, Phys. Rev. D71, 083502 (2005), astro-ph/0408426. 
[33] T. Kanzaki, M. Kawasaki, K. Kohri, and T. Moroi, Phys. Rev. D75, 025011 (2007), hep- 

34 



ph/0609246. 

[34] M. Kawasaki, K. Kohri, T. Moroi, and A. Yotsuyanagi, Phys. Rev. D78, 065011 (2008), 
0804.3745. 

[35] D. Cumberbatch et al., Phys. Rev. D76, 123005 (2007), 0708.0095. 
[36] M. H. Reno and D. Seckel, Phys. Rev. D37, 3441 (1988). 

[37] S. Dimopoulos, R. Esmailzadeh, L. J. Hall, and G. D. Starkman, Astrophys. J. 330, 545 

(1988) . 

[38] S. Dimopoulos, R. Esmailzadeh, L. J. Hall, and G. D. Starkman, Phys. Rev. Lett. 60, 7 (1988). 
[39] S. Dimopoulos, R. Esmailzadeh, L. J. Hall, and G. D. Starkman, Nucl. Phys. B311, 699 

(1989) . 

[40] M. Y. Khlopov, Y. L. Levitan, E. V. Sedelnikov, and I. M. Sobol, Phys. Atom. Nucl. 57, 1393 
(1994). 

[41] E. V. Sedelnikov, S. S. Filippov, and M. Y. Khlopov, Phys. Atom. Nucl. 58, 235 (1995). 
[42] K. Jedamzik, Phys. Rev. Lett. 84, 3248 (2000), astro-ph/9909445. 
[43] K. Jedamzik, Phys. Rev. D70, 063524 (2004), astro-ph/0402344. 
[44] K. Jedamzik, Phys. Rev. D70, 083510 (2004), astro-ph/0405583. 

[45] K. Jedamzik, K.-Y. Choi, L. Roszkowski, and R. Ruiz de Austri, JCAP 0607, 007 (2006), 
hep-ph/0512044. 

[46] K. Jedamzik, Phys. Rev. D74, 103509 (2006), hep-ph/0604251. 

[47] M. Kusakabe, T. Kajino, and G. J. Mathews, Phys. Rev. D74, 023526 (2006), astro- 
ph/0605255. 

[48] M. Kusakabe et al., Phys. Rev. D79, 123513 (2009), 0806.4040. 

[49] M. Pospelov and J. Pradler, Phys. Rev. D82, 103514 (2010), 1006.4172. 

[50] M. Pospelov and J. Pradler (2010), 1010.4079. 

[51] J. Melendez and I. Ramirez, Astrophys. J. 615, L33 (2004), astro-ph/0409383. 

[52] M. Asplund, D. L. Lambert, P. E. Nissen, F. Primas, and V. V. Smith, Astrophys. J. 644, 

229 (2006), astro-ph/0510636. 
[53] F. Spite and M. Spite, Astron. Astrophys. 115, 357 (1982). 

[54] S. G. Ryan, T. C. Beers, K. A. Olive, B. D. Fields, and J. E. Norris, Astrophys. J. 530, L57 

(2000), arXiv:astro-ph/9905211. 
[55] P. Bonifacio, P. Molaro, T. Sivarani, R. Cayrel, M. Spite, F. Spite, B. Plez, J. Andersen, 

35 



B. Barbuy, T. C. Beers, et al, Astron. Astrophys. 462, 851 (2007), arXiv:astro-ph/0610245. 

[56] J. R. Shi, T. Gehren, H. W. Zhang, J. L. Zeng, and G. Zhao, Astron. Astrophys. 465, 587 
(2007). 

[57] W. Aoki et al., Astrophys. J. 698, 1803 (2009), 0904.1448. 

[58] R. H. Cyburt, B. D. Fields, and K. A. Olive, JCAP 0811, 012 (2008), 0808.2818. 

[59] D. Larson et al. (2010), 1001.4635. 

[60] O. Richard, G. Michaud, and J. Richer, Astrophys. J. 619, 538 (2005), astro-ph/0409672. 

[61] A. J. Korn et al, Nature 442, 657 (2006), astro-ph/0608201. 

[62] K. Lind, F. Primas, C. Charbonnel, F. Grundahl, and M. Asplund (2009), 0906.2876. 

[63] R. Cayrel et al. (2007), 0708.3819. 

[64] M. Steffen, R. Cayrel, P. Bonifacio, H. G. Ludwig, and E. Caffau (2009), 0910.5917. 

[65] N. Prantzos, Astron. Astrophys. 448, 665 (2006), arXiv:astro-ph/0510122. 

[66] R. N. Boyd, C. R. Brune, G. M. Fuller, and C. J. Smith (2010), 1008.0848. 

[67] R. H. Cyburt and M. Pospelov (2009), 0906.4373. 

[68] N. Chakraborty, B. D. Fields, and K. A. Olive (2010), 1011.0722. 

[69] N. Arkani-Hamed and S. Dimopoulos, JHEP 06, 073 (2005), hep-th/0405159. 

[70] N. Arkani-Hamed, S. Dimopoulos, G. F. Giudice, and A. Romanino, Nucl. Phys. B709, 3 

(2005), hep-ph/0409232. 

[71] S. Raby, Phys. Lett. B422, 158 (1998), hep-ph/9712254. 

[72] S. Shirai, M. Yamazaki, and K. Yonekura, JHEP 06, 056 (2010), 1003.3155. 

[73] L. Covi, M. Olechowski, S. Pokorski, K. Turzynski, and J. D. Wehs (2010), 1009.3801. 

[74] U. Sarid and S. D. Thomas, Phys. Rev. Lett. 85, 1178 (2000), hep-ph/9909349. 

[75] J. Hisano, K. Nakayama, S. Sugiyama, T. Takesako, and M. Yamanaka, Phys. Lett. B691, 46 

(2010), 1003.3648. 

[76] K. Nakayama, F. Takahashi, and T. T. Yanagida (2010), 1010.5693. 

[77] J. Kang, M. A. Luty, and S. Nasri, JHEP 09, 086 (2008), hep-ph/0611322. 

[78] H. Baer, K.-m. Cheung, and J. F. Gunion, Phys. Rev. D59, 075002 (1999), hep-ph/9806361. 

[79] E. W. Kolb and M. S. Turner, The early universe (Addison- Wesley, 1990). 

[80] S. Wolfram, Phys. Lett. B82, 65 (1979). 

[81] C. B. Dover, T. K. Gaisser, and G. Steigman, Phys. Rev. Lett. 42, 1117 (1979). 

[82] G. D. Starkman, A. Gould, R. Esmailzadeh, and S. Dimopoulos, Phys. Rev. D41, 3594 (1990). 



36 



[83] K. Nakamura and Particle Data Group, Journal of Physics G Nuclear Physics 37, 075021 
(2010). 

[84] N. Prantzos, M. Casse, and E. Vangioni-Flam, Astrophys. J. 403, 630 (1993). 

[85] R. Ramaty, B. Kozlovsky, R. E. Lingenfelter, and H. Reeves, Astrophys. J. 488, 730 (1997), 

arXiv:astro-ph/9610255. 
[86] M. Kusakabe, Astrophys. J. 681, 18 (2008), 0803.3401. 
[87] S. E. Woosley and T. A. Weaver, Astrophys. J. Suppl. 101, 181 (1995). 

[88] T. Yoshida, T. Kajino, and D. H. Hartmann, Phys. Rev. Lett. 94, 231101 (2005), astro- 
ph/0505043. 

[89] E. Hiyama, M. Kamimura, T. Motoba, T. Yamada, and Y. Yamamoto, Phys. Rev. C66, 

024007 (2002), nucl-th/0204059. 
[90] E. Hiyama, Y. Kino, and M. Kamimura, Prog. Part. Nucl. Phys. 51, 223 (2003). 
[91] J. MartoreU, D. W. L. Sprung, and D. C. Zheng, Phys. Rev. C51, 1127 (1995). 
[92] A. Amroun et al., Nucl. Phys. A579, 596 (1994). 
[93] I. Tanihata et al., Phys. Lett. B206, 592 (1988). 

[94] M. Fukuda, M. Mihara, T. Fukao, S. Fukuda, M. Ishihara, S. Ito, T. Kobayashi, K. Matsuta, 

T. Minamisono, S. Momota, et al., Nuclear Physics A 656, 209 (1999). 
[95] C. Angulo et al, Nucl. Phys. A656, 3 (1999). 

[96] R. N. Boyd, An Introduction to Nuclear Astrophysics (the University of Chicago Press, 2008). 
[97] G. R. Caughlan and W. A. Fowler, Atom. Data Nucl. Data Tabl. 40, 283 (1988). 
[98] M. Sowerby, Journal of Nuclear Energy 24, 323 (1970). 

[99] C. A. Bertulani, Comput. Phys. Commun. 156, 123 (2003), nucl-th/0305099. 
[100] G. G. Simon, C. Schmitt, and V. H. Walther, Nucl. Phys. A364, 285 (1981). 
[101] D. R. Tilley et al, Nucl. Phys. A708, 3 (2002). 

[102] W. A. Fowler, G. R. Caughlan, and B. A. Zimmerman, Ann. Rev. Astron. Astrophys. 5, 525 
(1967). 

[103] L. Kawano, NASA STI/Recon Technical Report N 92, 25163 (1992). 
[104] M. S. Smith, L. H. Kawano, and R. A. Malaney, Astrophys. J. Suppl. 85, 219 (1993). 
[105] P. Descouvemont, A. Adahchour, C. Angulo, A. Coc, and E. Vangioni-Flam (2004), astro- 
ph/0407101. 

[106] G. J. Mathews, T. Kajino, and T. Shima, Phys. Rev. D71, 021302 (2005), astro-ph/0408523. 

37 



[107] Y. I. Izotov and T. X. Thuan, Astrophys. J. 710, L67 (2010), 1001.4440. 
[108] E. Aver, K. A. Olive, and E. D. Skillman, JCAP 1005, 003 (2010), 1001.5218. 
[109] M. Pettini, B. J. Zych, M. T. Murphy, A. Lewis, and C. C. Steidel (2008), 0805.0594. 
[110] T. M. Bania, R. T. Rood, and D. S. Balser, Nature 415, 54 (2002). 

[Ill] M. Kusakabe, T. Kajino, and G. J. Mathews, submitted to lOP Pubhsing to be pubUshed 
in INPC 2010 proceedings in the Open Access Journal of Physics: Conference Series. 

[112] M. Kamimura, Y. Kino, and E. Hiyama, in American Institute of Physics Conference Series, 
edited by H. Susa, M. Arnould, S. Gales, T. Motobayashi, C. Scheidenberger, & H. Utsunomiya 
(2010), vol. 1238 of American Institute of Physics Conference Series, pp. 139-144. 

[113] E. Rollinde, E. Vangioni, and K. A. Olive, Astrophys. J. 651, 658 (2006), astro-ph/0605633. 

[114] E. Rollinde, D. Maurin, E. Vangioni, K. A. Olive, and S. Inoue, Astrophys. J. 673, 676 
(2008), 0707.2086. 

[115] R. N. Boyd and T. Kajino, Astrophys. J. 336, L55 (1989). 

[116] T. Kajino and R. N. Boyd, Astrophys. J. 359, 267 (1990). 

[117] T. Kajino, G. J. Mathews, and G. M. Fuller, Astrophys. J. 364, 7 (1990). 

[118] A. Coc, P. Delb our go- Salvador, F. de Oliveira, P. Aguer, S. Barhoumi, G. Bogaert, J. Kiener, 
A. Lafebvre, and J. P. Thibaud, Astrophys. J. 402, 62 (1993). 

[119] M. Orito, T. Kajino, R. N. Boyd, and G. J. Mathews, Astrophys. J. 488, 515 (1997), 
arXiv:astro-ph/9609130. 

[120] C. Chiappini, A. Renda, and F. Matteucci, Astron. Astrophys. 395, 789 (2002), astro- 
ph/0209240. 

[121] E. Vangioni-Flam, K. A. Olive, B. D. Fields, and M. Casse, Astrophys. J. 585, 611 (2003), 

astro-ph/0207583. 
[122] K. Lodders, Astrophys. J. 591, 1220 (2003). 

[123] H. Ito, W. Aoki, S. Honda, and T. C. Beers, Astrophys. J. 698, L37 (2009), 0905.0950. 
[124] F. Hammache et al. (2010), 1011.6179. 

[125] M. Yahiro, Y. Iseri, H. Kameyama, M. Kamimura, and M. Kawai, Progress of Theoretical 

Physics Supplement 89, 32 (1986). 
[126] N. Austern et al., Phys. Rept. 154, 125 (1987). 

[127] Recently, a new measurement of the cross section of radiative a ca ptur e by deuteron and the 



-"Li abundance predicted based upon the result have been reported. 

38 



1241 ] 



[128] [A/H]=log(A/H)— log(A/H)0 is the number ratio of nuclide A to H measured in a logarithmic 

scale normalized to the solar value. 
[129] These parameters have been adjusted to fit the deuteron binding energy and low-energy 



125 



126 ] 



triplet-even proton-neutron scattering phase shifts 

[130] This value is an example which leads to the binding energy of deuteron when another pa- 
rameter is fixed to be fow = —25.5 MeV. 

[131] The ^Bex nucleus produced in this reaction pathway immediately p-decays into ^Hex and 
two protons. 

[132] This reaction can be efficient since the reaction Q-value of X(Jhi, t)^Hex is very small (see 
Table HD. 



[133] WWW: http://www.tunl.duke.edu/nucldata/index.shtml 



[134] The reaction rates and related information of ^llex{d,^hi)X , ^Hev (t7Li)X and 



^Hejis:(^He/Be)X in the case of leptonic X particle can be found in Refs. [8 



15| 



39 



